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We investigate the combined influence of a magnetic field and a harmonic interparticle interac- 
tion on the thermodynamic properties of a finite number of spin polarized fermions in a confiment 
potential. This study is an extension using our path integral approach of symmetrized density ma- 
, trices for identical particles. The thermodynamical properties are calculated for a three dimensional 

model of N harmonically interacting spin polarized fermions in a parabolic potential well in the 
presence of a magnetic field. The free energy and the internal energy are obtained for a limited 
number of particles. Deviations from the thermodynamical limit become negligible for about 100 
or more particles, but even for a smaller number of fermions present in the well, scaling relations 
similar to those of the continuum approximation to the density of states are already satisfied. 



PACS: 05.30.-d, 03.75.Fi, 32.80.Pj. 



I. INTRODUCTION 



C 

In the present paper we study the thermodynamical properties of a confined system of spin-polarized fermions in 
i ^ i ■ the presence of a magnetic field. The method used is an extension of the combination of the path integral formalism 
jjj and the method of symmetrized density matrices ^ , developed previously for a model system of harmonically 
t-H , interacting identical particles (bosons or fermions) in a parabolic well (hereafter for brevity referred to as the harmonic 
^ model). 

Because of the experimental realization of Bose-Einstein condensation |@-^] and the theoretical work on this phe- 
nomenon employing other methods |10-17j, full details for this harmonic model with interaction were first worked out 



for bosons. The model shows the onset of Bose-Einstein condensation in the specific heat |18[ for a finite number of 
particles, and its moment of inertia is drastically reduced below the condensation temperature JTsj ] . An application 
Q"\ ' of the method to real systems can be found in Ref. §2j§ for 87 Kb. 

0^ ■ The actual calculations for the fermion case require more advanced techniques, such as the generating function 
approach and the corresponding contour integration, because of a numerical sign problem. In the absence of a 
magnetic field, explicit results for the thermodynamics and the static correlation functions of the harmonic model of 
spin-polarized fermions were already obtained with these techniques |^|| . 

The harmonic model clearly has intrinsic value on its own, because it is one of the rare examples of an exactly 
soluble many-particle system with interaction. The physics of the model is relatively straightforward in the sense 
that it allows for center of mass excitations that oscillate at frequencies different from those of the internal degrees of 
freedom. This property makes it well suited as a trial model for the variational treatment of the thermodynamics of 
systems with more realistic interactions because the model parameters can be related with the system characteristics 
with the aid of the Jensen- Feynman inequality . The present paper adresses only the first part of such an approach 
' because it requires also the density and the pair correlation function which we could obtain under the simplifying 
assumptions of no magnetic field Furthermore, it provides a testing ground for new approaches to Monte Carlo 
simulations of interacting fermions such as many body diffusion pfi-G3|. Especially for quantum dots, it is important 



to take the magnetic field into account in order to freeze out the opposite spin states. In the present paper we present 
an extension of the methods mentioned above to harmonically interacting confined fermions in a magnetic field. 



*Also at Universiteit Antwerpen (RUCA), and at Technische Universiteit Eindhoven, NL 5600 MB Eindhoven, The 
Netherlands. 
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The paper is organized as follows. In Sec. II we present the path integral for harmonically interacting particles in a 
parabolic confinement potential in the presence of a homogeneous magnetic field. This will be done for distinguishable 
as well as for identical particles. The mathematical details of the calculation for identical particles are given in 
Appendix A. In Sec. Ill the permutation symmetry will be taken into account with the aid of the projection technique. 
The introduction of the permutation symmetry implies the rewriting of the sum over all possible permutations to a 
cyclic summation Q which leads to the generating function of the partition function. Specific results for fermions will 
be presented in Sec. IV. This involves the extraction of the partition function and other thermodynamical quantities 
from the generating function. Also the ground state energy and the magnetic susceptibility in the zero-temperature 
limit will be investigated. Additionally we will study the finite number corrections to the thermodynamic limit for 
the free energy and the internal energy as a function of temperature and magnetic field. In the last section some 
conclusions are given. 



II. IDENTICAL OSCILLATORS IN A MAGNETIC FIELD 



The calculation of the path integral for N identical interacting oscillators in a magnetic field is similar to the case 
without a magnetic field in || . This approach crucially relies on the detailed investigation of the classical action, 
and the path-integral corrections to this classical action. Alternatively, a more stochastic approach could be followed 
p4| . The Lagrangian (in atomic units) for N oscillators with harmonic two-body interactions and in the presence of 
a homogeneous magnetic field is given by 

L = -^(rj-2^)-^-^ (1) 

=1 



where uj c is the cyclotron frequency and 

O 2 N ,? N 

v, = \Y.^ and ^±TE( r rn) 2 (2) 

3=1 J,2=l 

It is obvious that the two-body potential is either attractive or repulsive depending on the plus sign or the minus sign 
considered in Vi. The magnetic field introduces a coupling in the plain perpendicular to its direction. This means 
that we can separate the Lagrangian into two contributions L = L xy + L z . The Lagrangian L z simply describes a 
harmonic oscillator, whereas L xy contains the magnetic field. The Lagrangian can be rewritten in terms of the center 
of mass coordinate R (X, Y, Z) and the coordinates r/j (uj , Vj , Wj ) describing the coordinates of the particles measured 
from the center of mass 

1 N 

R =T^ r J' %= r i- R > (3) 



from which 



N ■ 



i 2 N 

v 1 + v 2 = v CM + v, v C M = 2 Nn2R2 > v = T^A 2 ' (4) 

3=1 



with 



w 



yJVL 2 ± Nuj 2 . (5) 



For a repulsive two-particle potential the internal frequency w has to satisfy the stability condition that the confining 
potential has to be sufficiently strong to overcome the repulsion between the particles. We draw attention to the 
fact that the transformation to the center of mass coordinate system diagonalizes neither the Lagrangian nor the 
Hamiltonian, because of the subsidiary condition 



N 



5> 3 --R) = 0. (6) 
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We obtain the propagator for distinguishable (indicated by a subscript D) particles from the action expressed in 
the imaginary time variable = X/kT and it is written as 

K D (r'/, . . . , r£, | ri, . . . , r' N , 0) = K D ((x", y") , | (x', y') , 0) x Kjj (z", \ z' , 0) , (7) 

where the vector x denotes the N dimensional x coordinates of the particles, with the notation x T = (x\, x%, . . . , Xn) , 
and similarly for y and z. The propagator for a single oscillator with frequency w in a magnetic field is well known 
H and given by 



ifW(r",/3|r',0) =J 



2tt sinh 0w 2tt sinh 0s 



x exp ■ 



x exp ■ 



(V) 2 + (z") 2 ) cosh/3' 



w - 2z z 



2 s'mh0za 

s ([x"f + [y"f + {x'f + (y') 2 ) cosh /3s - 2 (x'x" + y'y") cosh/3w L 
2 sinh /3s 



x rx], {-,[ co L (x"y" ~ x'y') - s*-^^ (y'x" ~ *V) ) ^ . (si 



where u/, = ^ is the Larmor frequency and the eigenfrequency s is given by 



2+ w £. (9) 

The propagator for TV distinguishable interacting oscillators in a magnetic field thus becomes 

K n (VNZ",(3\ VNZ',o) K ULjSCM (VnX",VnY",0\ ^NX',VNY',o) 
K d (r", (3 |f',0) ) f -A 

k w (y/NZ", I Vnz', oj /y ui , s (Vivx", VivY", I VFa-', VFy , oj 

AT 

x H K UL , S (x'!, y'J, | xj, yj, 0) K w (z'J, \ , (10) 

3=1 

where 



s = y w 2 + w| and s C m = y ^ 2 + w£. (11) 

The factor \/]V in the center-of-mass coordinates in (|l0|) describes the mass iV (in atomic units) of the center. The 
denominator in ( |Io| ) accounts for the fact that the internal degrees of freedom are linearly dependent because of the 
subsidiary conditions (|^). Intuitively this factor is quite natural, because the propagator would be the product of 
one-particle propagators || if the particles were independent. 

Knowing the propagator for distinguishable particles, the symmetrized density matrix Ki for identical particles can 
be obtained through the appropriate symmetric or antisymmetric projection 

p 

where P denotes the permutation matrix, with £ = +1 for bosons and £ = — 1 for fermions. Even for this harmonic 
model, with or without magnetic field, the sum over the permutations has to remain rather formal at the level of the 
propagator. However, for the partition function 

Zi (0, N) = J dm (f,0 I r, 0) = i Z P J dfR D /? I r, 0) , (13) 

analytical progress can be made with this summation as will be discussed in the next section. First of all, one has to 
deal with the center-of-mass contribution to the propagator. Afterwards, the summation over all possible permutations 
will be rewritten as a summation over all possible cycles. 
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III. GENERATING FUNCTION OF THE PARTITION FUNCTION 



The center of mass is not independent of the positions of the other particles, which complicates the calculation of 
the trace of the propagator. To deal with the contribution of the center-of-mass coordinate R to the propagator, we 
introduce the delta function S (li — J2f=i r i) m its Fourier representation as in ||. This delta function allows one 
to formally treat the center-of-mass coordinate as an independent variable. Applying this identity to the partition 
function, one ends up with 

r r dk , D K a {^z, 1 Vnz, o) k ul scm (Vnx, Vny, p \ Vnx, Vny, o) 

Zi((3,N) = dR --^e lkR ) f : -A ^ 

J J (2ir) K... I \f~N 7, 8 I \fW ZJn K I \/N X \/NY. 3 I \f~N X \/~NY. () I 



( 2 ^) 3 K w (y/NZ, 13 | VNZ, o) K^ l S (VNX, VNY, (3 | y/NX, VNY, o) 

r i N 

x / dT m E f II k "l,> (( Px h • ( p yh > i x i> %>°) K -> {( pz h >0 l °) e- i% - p *' N . (w) 
J ■ p j=x 

The problem at hand is the rewriting of the summation over the permutations as a sum over all possible cycles [pj. 
This cyclic decomposition requires the solution of the path integral for a driven harmonic oscillator in a magnetic 
field, which is discussed in appendix A. 

A permutation can be decomposed into Mi cycles of length £, and the positive integers Mg and I have to satisfy 
the constraint 

^tM t = N. (15) 

i 

The number M (Mi, . . . ,Mn) of permutations with Mi cycles of length 1,. . ., Mg cycles of length I,... is given 
by M (Mi, . . . , Mat) = Nl/ \J\ e Mgl£ Ml ] . Furthermore, a cycle of length I will be obtained from t — 1 permutations. 
Thus, the sign factor £ p can be rewritten as £ p = Y[g £^ _1 ) Mf . These considerations enable one to rewrite the partition 
function as 

r f dk D K a (Vnz, P I Va?z, o) k ul scm (Vnx, Vny, 13 I Vnx, Vny, o) 

Zj(/3,N) = dR ^V kR ) L ■ ^ J- 

■' J ( 27r ) K w (VN Z, (3 | VN Z, OJ K ULiS (yNX,VNY,p\VNX,VNY,0\ 

M e \e M - 



N M-l)M e 

E n IZITmT [M k )r < (16) 

Mi,...,Mn j=l 



with 



K t (k) = fdvt+i... fdviS(r £+ i -v^flK^s (x j+ i,y j+ i,(3 \ Xj ,y 3 ,Q)K w (z J+1 ,(3\ Zj ,0) e -^/ N . (17) 
J J j=i 

The delta function explicitly indicates that the trace is taken over a cycle of length I. It is obvious that /Q (k) 
factorizes as 

Kt(k)=K, i (k x ,k v )K t (k z ). (18) 

Taking into account the semigroup property of the propagators K ULtS (xj+i,yj+i, (3 \ Xj,yj, 0) and K w (zj+i,{3 \ Zj, 0), 
one immediately recognizes in Kg. (k) the partition function of a driven harmonic oscillator in a magnetic field, with 
the driving force 

*(T)=ijjY,&<J-3p). (19) 
j=o 

The calculation of the propagator and the partition function for a driven harmonic oscillator in a magnetic field, with 
the Lagrangian 

4% = \i ±2 + i 2 ) - 2uJ ^xy - ^ (x 2 + y 2 ) + f x (t) x + f v (r) y, (20) 
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is illustrated in appendix A. As mentioned above, the hard core of this approach is the evaluation of the classical 
action, but a fully stochastic method (24} could as well have been followed. After tedious algebra one eventually finds 
for the partition function 

, K n (Vnz, p I Vnz, o) k ul scm (Vnx, Vny, p I Vnx, Vny, 

Z T {P,N) = / dK — ) (■ — : — J> r 

J k w [Vnz, p | Vnz, oj k ULi . [Vnx, Vny, p \ Vnx, Vny, oj 

dk ik-R e ( 1 k l + k l sinh /3s 1 k\ sinh/?w 



(27r) \ 4s N cosh /3s — cosh /3w_l 4w N cosh f3w — 1 / 

X M^M* ,9 ( 8 Sinh ^ Sinh Sinh ¥ J ' (21) 

The remaining integrations over k and R are Gaussian and relatively easy to perform, leading to 

sinh r ^pl sinh sinh 2? 

with 

Z/ N) = M^M, 5 ' 3 ^ ]Sr V 8 sinh ^±^1 sinh sinh ^ J ' ^ 

The contribution Zj (/3, N) derives from the internal degrees of freedom, treated as independent particles. It contains 
the full influence of the statistics of the particles, and leads to the true partition function Zi (/3, TV) by multiplication 
with a simple analytical factor. In practice, the condition ( |l5| ) complicates the use of the above expression for the 
partition function for a large number of particles. However this difficulty can be overcome through the use of the 
generating function. From the generating function one can then extract the partition function through an inversion 
of its defining Taylor series. 



A. Generating function and recurrence relation for the partition function 

The generating function technique was used before || to obtain the partition function of a set of harmonically 
interacting identical oscillators in the absence of a magnetic field. In the presence of a magnetic field, a similar 
construction can be used. Introducing the generating function as 

oo 

3 J (l3,u)=Y i Z I {N)u N , (24) 

JV=0 

(with Zj (/3, 0) = 1 by definition), the partition function for the internal degrees of freedom can be obtained from 

1 d N 

z ' ( ^ ) = m^ SjWI «=°- (25) 

The generating function itself can be obtained with straightforward algebra 

M ft„)=exp(|:— _ I __ n __ I j, (26) 

with the following notation 

b = e- /3w , b 1 = e- fj{s+UJL \ b 2 = e-^ s - UL \ (27) 

The cyclic summation can be rewritten in terms of the occupation number representation which directly involves the 
single particle energy levels: 
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(28) 



By simply applying the chain rule, the expression (25) for the partition function can be written as a recurrence relation 



JV 



Zj (/3,7V) 



(1-6* 



(bib 2 by 



b{) (1-60(1 



-Z/ (/3, N -I). 



(29) 



However if the number of particles increases, this recurrence relation becomes numerically unpractical because of a 
numerical sign problem for fermions and drastically increasing simulation time for bosons. For the remaining part of 
this paper, the attention will be focussed on the fermion case. 



IV. THERMODYNAMIC AL PROPERTIES 



The thermodynamical properties of the fermion model can in essence be determined from the contribution Zp ((3, N) 
of the internal degrees of freedom. As is clear from (^2|), the center of mass correction only adds a trivial contribution to 
the free energy. All the effects of the fermion statistics are collected in Zp (f3, N) . We first study the zero temperature 
limit, in which special attention will be paid to the two dimensional case in the xy-plane, and subsequently the 
evolution of the free energy and the internal energy as a function of the temperature and of the magnetic field. 



A. Zero temperature limit 



The ground state properties of the fermion model crucially depend on the single-particle energy levels E VMl ^ 2 . 
These levels and their occupation by fermions have been discussed and plotted earlier, e.g., in Ref. p5j. For easier 
reference, we plot the 20 lowest levels Fig. 1. To guide the eye, the Fermi energies Ep corresponding to the fully 
occupied levels at ujl = arc indicated by the dashed line. Note that the magnetic field does not substantially 
influence the magnitude of the Fermi energy, which remains of order TV 1 / 3 for sufficiently large N. The magnetic 
field immediately lifts the degeneracy, but with increasing magnetic field other degeneracies appear and disappear 
again at particular values of the magnetic field. Although these degeneracies have little effect on the magnitude of the 
ground state energy Eg = ^2e<e f ^, they have a drastic effect on the magnetic susceptibility, which is proportional 
to dEa/dujL, as shown in Fig. 2 as a function of the magnetic field. The discontinuities in the magnetic susceptibility 
occur at those values of the magnetic field where the single-particle energies become degenerate. 



B. Free energy and internal energy 



As mentioned above, the sign problem for fermions can be worked around by inverting the defining Taylor-series 
( pif ) for the generating function. The Fowler-Darwin method [g6j provides an accurate and elegant way to realize 
this inversion: 

z,(M0 = ^jf ^^d*. (30) 

If one considers a circular contour z — ue l ° with radius u, an optimal value of u can be determined by the method of 
steepest descent 

— (In Sj? (13, u) - Ainu) = N = u-^-lnZp ((3,u). (31) 
du du 

Using (p8|), this condition becomes 

oo 

N= Yl n»w (32) 

with 



G 



1 + e 



0E„ 



= VI 



si 



S2 



(33) 



If u would be interpreted as the fugacity u — e^** with chemical potential fi, one thus would recover similar results 
as for the expectation value of the number of particles in the grand canonical ensemble. The result for the chemical 
potential as a function of temperature is shown in Fig. 3. for various values of the magnetic field and for N = 2. 
The chemical potential is plotted in units of the chemical potential at T = 0, and the temperature in units of wN 1 ^ 3 , 
which is the order of magnitude for the Fermi energy. In Fig. 4 and 5 the corresponding results are shown for N = 10 
and N — 100. For N > 100 it turns out that /z (T) j \i (T = 0) as a function of kT/wN 1 / 3 becomes almost independent 
of both the number of particles and of the magnetic field. 

However, in the present treatment the determination of u = from ( |32| ) only gives the zero-order contribution 
to the partition function. A correction by the integration factor in ( |30| ) has to be applied. Using the symmetry of the 
integrand in (|30|), the partition function can be rewritten as 



If (13, N) 



E F (J3,u) 1 f 27r 5f (P,ue ie ) 
Z F ((3,u) 



-iN6 



with 



2vr 



^ (6») = Re 



dO = 



~ F (f3,u) 



,N 



* (0) dO, 



1 



-i8N 



5f (A 



3 f {fi, u) 



(34) 



(35) 



The function if? (9) has to be calculated and integrated numerically. The determination of the free energy 

r(°) <i\ x i _ .". i,, / \T; i(4\ ,m (35) 



F F (/3, N) = ~ lnZ F (fi, N) = ((3, N) - \ In / * (0) d0 
^ P Jo 



then becomes straightforward, with the zero-order contribution F^P (/?, AT) 



fl 111 U N 



from the steepest-descent 
For = 2w the free energy 



approximation. In the absence of a magnetic field, the results are discussed in |^ 
per particle in units of the fermi energy is plotted as a function of kT/E F for 10 fermions in Fig. 6 and for 100 
fermions in Fig. 7, and compared to the zero-order steepest descent contribution. Again, for N > 100 the finite 
number corrections upon the thermodynamical limit become negligible for all practical purposes. The internal energy 
Up = -jp (fi¥p) shows the same universality, as is shown in Fig. 8 where the internal energy per particle in units of 
the Fermi energy is plotted versus the temperature in units of the Fermi temperature for lul — 0, w and 2w. 



V. CONCLUSION AND DISCUSSION 



Using the path integral approach of symmetrized density matrices for identical particles, the thermodynamical 
properties were calculated for a three dimensional model of N harmonically interacting spin-polarized fermions in a 
parabolic potential well in the presence of a magnetic field. The method used is a generalization of the procedure 
developed earlier in the absence of a magnetic field. Explicit results were obtained for the ground state energy, the 
free energy and the internal energy for a limited number of particles. The model can be described as a number of spin 
polarized identical particles in a parabolic confinement potential interacting through a special many body interaction 
with the consequence that the center of mass is allowed to move independent from the other degrees of freedom. For 
an analoguous model other forms of confinement potentials have been investigated without two body interaction. p9| ] 

The statistics with a finite number of particles in the confinement potential and the cross over to density dependent 
expressions known from the thermodynamical limit can be studied in this model: as soon as the number of fermions 
is sufficiently large (in the order of N > 100) the results are shown to agree and the finite number corrections become 
relatively small. The internal energy turns out to obey a scaling law, similar to the scaling from the continuum 
approximation for the density of states. 
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APPENDIX A: THE PATH INTEGRAL OF THE MODEL FOR DISTINGUISHABLE PARTICLES IN 
THE PRESENCE OF A MAGNETIC FIELD, AND A TIME-DEPENDENT DRIVING FORCE. 

The propagator of a two dimensional harmonic oscillator in the presence of a magnetic field, characterized by the 
Larmor frequency ul = w c /2, and under the influence of a time-dependent driving force f = {f x ,fy), provides the 
basic building blocks for the harmonic model system of identical interacting particles which is the subject of the 
present paper. Although the calculation of this propagator relies on standard techniques, to the best of our knowledge 
it is not documented in the literature. Therefore we discuss its derivation here in some detail. The Lagrangian under 
consideration is given by (in atomic units h = m = \e\ = 1) 

41 = \ + V 2 ) ^Lxy ^ (x 2 + y 2 ) + f x (r) x + f y (r) y. (Al) 

The classical equations of motion in Euclidean time r = it are 

d^ x du 

- — = -2iu; L -£-w 2 x + f x , (A2) 
d^u dx 

-^ = ^ L --w 2 y + f y . (A3) 

This set of coupled differential equations can be solved and yields 

x (t) = x h (t) + x p (t) , y (t) = y h (t) + y p (r) . (A4) 

The solutions x^ (r) and yh (r) of the homogeneous equations of motion (without the driving force) , which exhaust 
the boundary conditions Xh (0) = x' , Xh (0) = x", yh (0) = y' , yh (0) = y", are found to be 

x h (t) \ _ sinh,s(/3 — t) / cosh lolt isinhujLT \ ( x 
yh{r) J ~~ sinh/3s ^-isinh^LT cosher J \ y' 

sinh st / coshwi (0 - r) —i sinh uj l (0 - r) 
sinh 0s \ i sinh u) L {0-t) cosh uj l (0-t) 



(A5) 



with 



= \A 



u 2 L +w 2 . (A6) 



The derivation of the particular solutions x p (r) and y p (r) , with the boundary conditions x p (0) = x p (0) = y p (0) 
Up = 0, is slightly more involved but eventually results in 



x p (t) \ _ 1 sinh s (0 - r) ( T / coshw L (r-cr) i sinh lu l (r - a) \ ( f x (a) 
y P ( T ) J _ s sinh 0s J \ -i sinh u L (r - a) coshw L (r-cr) J \ f y (a) 



sinh scrda 



+ lsinW ffi/ cosh^^-a) (r - a) \ ff f)\ ^ a 

ssmh0sJ T \ -i sinh u L (r - a) coshw L (r-cr) / y f y (a) J 

Given the classical trajectory with initial position (x',y') and final position (x",y") after an imaginary time lapse 0, 
the corresponding classical action St tC i = J Q Lfdr can be found by elementary methods. This eventually results in 



S f , cl (x",y",0\x',y',O) 



s 



| ({x") 2 + (y" f + (x' f + (y') 2 ) cosh/3s 



sinh 0s \ - (x'x" + y'y") cosh/3w L + i {x'y" - x"y') smh0uj L 



iuj L (x"y" - x'y') 
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+ 



sinh (is 



( +x' Jq (fx (t) cosher — if y (r) sinhu^r) sinhs ((3 — r) dr \ 
+ x " Jg (fx ( T ) coshuj L (f3 - t) + ify (r) sinhw L ((3 - r)) sinh srdr 
+V' L (fy ( r ) coshw L r + if x (t) sinhu L r) sinhs (/3 - r) 

V +y" In (fv ( r ) coshw i - r) - i/ x (r) sinhw L (/? - r)) sinh srdr / 



1 

s sinh /3s 7 



18 /-r 



(/x (t) /a, (ct) + /j, (r) /j, (a)) coshw L (r - cr) 
+i (fx (t) f y (a) + f y (r) f x (a)) sinhw L (r - cr) 



sinh as sinh s(fl — t) dadr. 



(A8) 



Since the Lagrangian is quadratic in the coordinates and the velocities, the quantum mechanical propagator is 
determined by the classical action, apart from a trivial normalization factor. The latter can be determined by 
elementary methods. The calculation presents no difficulties and results in 



K f (x",y",p\x',y',0) 



2n sinh /3s 



exp [S f , c i(x",y",p\x',y',0)], 



(A9) 



If one takes the limit u>l — * 0, the correct result M is recovered. 

For the treatment of the cyclic summations for identical particles, it is essential to know the trace of this propagator 



Kf (x,y,/3\ x,y,0) dxdy. 



The calculation of this quantity is straightforward and after some algebra one obtains 



Z f 



1 



2 (cosh /3s — cosh (3ujl) 



exp 



4s (cosh /3s — cos1i/3u;l) 



(A10) 



(All) 



with 



$09) 



r P 



I I it(\fi\ t tt\t(\\( coshwL (r - cr) sinhs (/? — |r - a\) 

< o j Q (f x (T)f x (<r)+fy(T)f»(<r))i +C0Shw ;^_( T _ CT | ) v siml ^ (r _ 1 ;) 



/3 /-/3 



-'o 



(fAT)fy(<j) + fy(T)fA°)) 



sinh W£ (r — cr) sinh s([3 — |<7 — r|) 
— sinhcji (J3 —\t — cr|) sinh s (r — cr 



dadr 
dadr. 



(A12) 



Again, if one takes the limit of a vanishing magnetic field one finds the correct result 
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Figure captions 

Fig. 1: Lowest single particle energy levels (in units of w) as a function of the Larmor frequency. The Fermi 
energies corresponding to 1, 4, 10 and 20 fcrmions (i.e. for closed shells in the absence of a magnetic field) are 
emphasized by dashed lines. The results can also be found in Ref. 

Fig. 2: Scaled magnetic susceptibility (1/-/V) (dEc / dui l) as a function of the magnetic field for 1, 4, 10, 20 fermions 
in the ground state. 

Fig. 3: Scaled chemical potential /i (T) //i(T = 0) as a function of the scaled temperature t — kT/wN 1 / 3 for 2 
fermions and for lul/w — 0, 1 and 2. 

Fig. 4: Same as Fig. 4, but for 10 fermions. 

Fig. 5: Same as Fig. 4, but for 100 fermions. 

Fig. 6: Scaled free energy per particle / = Fp/NEp as a function of the scaled temperature kT/Ep for 10 fermions 
and with the Larmor frequency = 2w. The zero-order "steepest descent" contribution is indicated by the 
dashed line. 

Fig. 7 : Same as Fig. 6, but for 100 fermions. 

Fig. 8: Scaled internal energy per particle u = Vp/NEp for 100 fermions as a function of the scaled temperature 
kT / Ep for several values of the Larmor frequency u>l = 0, w and 2w. 
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Fig. 7 




Fig. 8 



